Optimal Control for Generating Quantum Gates in Open Dissipative Systems 
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Optimal control methods for implementing quantum modules with least amount of relaxative loss 
are devised to give best approximations to unitary gates under relaxation. The potential gain by 
optimal control using relaxation parameters against time-optimal control is explored and exemplified 
in numerical and in algebraic terms: it is the method of choice to govern quantum systems within 
subspaces of weak relaxation whenever the drift Hamiltonian would otherwise drive the system 
through fast decaying modes. In a standard model system generalising decoherence-free subspaces 
to more realistic scenarios, openGRAPE-derived controls realise a cnot with fidelities beyond 95% 
instead of at most 15% for a standard Trotter expansion. As additional benefit it requires control 
fields orders of magnitude lower than the bang-bang decouplings in the latter. 

PACS numbers: 03.67.-a, 03.67.Lx, 03.65.Yz, 03.67.Pp; 82.56.Jn 



I. INTRODUCTION 

Using experimentally controllable quantum systems to 
perform computational tasks or to simulate other quan- 
tum systems p], [3] is promising: by exploiting quan- 
tum coherences, the complexity of a problem may re- 
duce when changing the setting from classical to quan- 
tum. Protecting quantum systems against relaxation is 
therefore tantamount to using coherent superpositions as 
a resource. To this end, decoherence-free subspaces have 
been applied 0], bang-bang controls 0| have been used 
for decoupling the system from dissipative interaction 
with the environment, while a quantum Zeno approach 
[j| may be taken to projectively keep the system within 
the desired subspace @. Controlling relaxation is both 
important and demanding [7|, |8|, Is L 1 1(1 . also in view of 
fault-tolerant quantum computing [llj or dynamic error 
correction [13]. Implementing quantum gates or quan- 
tum modules experimentally is in fact a challenge: one 
has to fight relaxation while simultaneously steering the 
quantum system with all its basis states into a linear im- 
age of maximal overlap with the target gate. — Recently 
we showed how near time-optimal control by GRAPE [l3| 
take pioneering realisations from their fidelity-limit to the 
decoherence- limit [13 ] . 

In spectroscopy, optimal control helps to keep the 
state in slowly relaxing modes of the Liouville space 
[Till . fl7| . In quantum computing, however, the en- 
tire basis has to be transformed. For generic relaxation 
scenarios, this precludes simple adaptation to the entire 
Liouville space: the gain of going along protected dimen- 
sions is outweighed by losses in the orthocomplement. 
Yet embedding logical qubits as decoherence-protected 
subsystem into a larger Liouville space of the encoding 
physical system raises questions: is the target module 
reachable within the protected subspace by admissible 
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controls? 

In this category of setting, the extended gradient algo- 
rithm openGRAPE turns out to be particularly powerful 
to give best approximations to unitary target gates in 
relaxative quantum systems thus extending the toolbox 
of quantum control, see e.g. 0, [H, [H, US HH, H3, 0, 
HH ]25|- [2^1 . Moreover, building upon a precursor of this 
work [22}, it has been shown in j28(] that non-Markovian 
relaxation models can be treated likewise, provided there 
is a finite-dimensional embedding such that the embed- 
ded system itself ultimately interacts with the environ- 
ment in a Markovian way. Time dependent T(t) have 
recently also been treated in the Markovian [29I [30| and 
non-Markovian regime [3lj . 

Here we study model systems that are fully controllable 
[S3, US, [H,[35|], i.e. those in which — neglecting relaxation 
for the moment — to any initial density operator p, the en- 
tire unitary orbit U(p) := {V pU~ x \ U unitary} can be 
reached [36j by evolutions under the system Hamiltonian 
(drift) and the experimentally admissible controls. More- 
over, certain tasks can be performed within a subspace, 
e.g. a subspace protected totally or partially against re- 
laxation explicitly given in the equation of motion. 



II. THEORY 

Unitary modules for quantum computation require 
synthesising a simultaneous linear image of all the basis 
states spanning the Hilbert space or subspace on which 
the gates shall act. It thus generalises the spectroscopic 
task to transfer the state of a system from a given initial 
one into maximal overlap with a desired target state. 



A. Preliminaries 

The control problem of maximising this overlap sub- 
ject to the dynamics being governed by an equation of 
motion may be addressed by our algorithm GRAPE [13 ]. 
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For state-to-state transfer in spectroscopy, one simply 
refers to the Hamiltonian equations of motion known as 
Schrodinger's equation (for pure states of closed systems 
represented in Hilbert space) or to Liouville's equation 
(for density operators in Liouville space) 



P) = -iH\i>) 
p = -i[H,p] 



(1) 
(2) 



In quantum computation, however, the above have to 
be lifted to the corresponding operator equations, which 
is facilitated using the notations Adu(-) ■— U(-)W and 
ad H (•) := [H, (•)] with U := e~ ltH obeying 



-it ad 



H (0=Ad u (.) 
and using 'o ! for the composition of maps in 



U = 



-iH U 

—i adn ° Adjj 



(3) 



(4) 
(5) 



These operator equations of motion occur in two scenar- 
ios for realising quantum gates or modules U(T) with 
maximum trace fidelities: The normalised quality func- 
tion (setting N := 2" for an n-qubit system henceforth) 



/' := i Rctr{L/ t t argct (7(r)} 



(6) 



covers the case where overall global phases shall be re- 
spected, whereas if a global phase is immaterial (22| 
(while the fixed phase relation between the matrix 
columns is kept as opposed to ref. [S3|), the quality func- 
tion 



/:=iRetr{Ad^ 



AH 



}=|/'l 



(7) 



applies. The latter identity is most easily seen [22j in 
the so-called wee-representation of p where one gets 
the conjugation superoperator Adjj = U (g> U (with U 
denoting the complex conjugate) and the commutator 
superoperator ad# = 1 ® H — H l <g> 1. 



B. Open GRAPE 

Likewise, under relaxation introduced by the operator 
r (which may, e.g., take GKS-Lindblad form), the re- 
spective Master equations for state transfer [35|] and its 
lift for gate synthesis read 



p = -(iadjj +T) p 
F = -(iad H +T) o F 



(8) 
(9) 



Again with N := 2" in n-qubit system, F denotes a quan- 
tum map in GL(N 2 ) as linear image over all basis states 
of the Liouville space representing the open system. The 
Lie-semigroup properties of F(t) have recently been elu- 
cidated in detail [33|: it is important to note that only in 
the special (and highly unusual) case of [&dn , F ] = the 



po = Pse(0) ® Pb(0) 
pse(O) 



AdwW > p{t) = W{t)p W\t) 
n SB I tr B 

ps sit) 
n s |tr E 
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FseW 



F S (t) 



ps(0) 



Figure 1: Time-evolution of a quantum system (S) embed- 
ded in some environment (E) and coupled to a bath (B). If 
the universal system evolves under the global unitary W(t), 
openGRAPE provides optimal controls for the category of sce- 
narios for which there is a finite-dimensional embedding such 
that the embedded system (SE) follows a time evolution un- 
der a Markovian quantum map FsE(t), while the reduced sys- 
tem of concern (S) may evolve in a non-Markovian way by a 
generic quantum map Fs{t). In the simplest of cases, Psb(O) 
is of tensor product form with Fs(t) being Markovian itself. 



map F(t) boils down to a mere contraction of the unitary 
conjugation Ady. In general, however, one is faced with 
an intricate interplay of the respective coherent (zad#) 
and incoherent (r) part of the time evolution: it explores 
a much richer set of quantum maps than contractions of 
Ad;/, as expressed in [39( in terms of a I, p-decomposition 
of the generators in g[(iV 2 ,C) of quantum maps. As will 
be shown below, it is this interplay that ultimately en- 
tails the need for relaxation-optimised control based on 
the full knowledge of the Master Eqn. Q, while in the 
special case of mere contractions of Adu, tracking maxi- 
mum qualities against fixed final times ('top curves', vide 
infra, e.g. Fig. [3] (a) upper panel) obtained for T = plus 
an estimate on the eigenvalues of T suffice to come up 
with good guesses of controls. 

Now for a Markovian Master equation to make sense in 
terms of physics, it is important that the quantum sub- 
system of concern is itself coupled to its environment in a 
way justifying to neglect any memory effects. This means 
the characteristic time scales under which the environ- 
ment correlation functions decay have to be sufficiently 
smaller than the time scale for the quantum evolution 
of the subsystem (see, e.g., [ifl|). — More precisely, as 
exemplified in Fig.[TJ we will assume that either the quan- 
tum system itself (S) or a finite-dimensional embedding 
of the system (SE) can be separated from the environ- 
mental bath (B) such that (at least) one of the quantum 
maps of the reduced system Fs(t) or FsE(t) is Marko- 
vian and allows for a description by a completely posi- 
tive semigroup (ill . l42l |43| , if the time evolution for the 
universal composite of (embedded) system plus bath is 
unitary. Examples where Fs(t) is Markovian have been 
given in a precursor (22j to this study, while a concrete 
setting of a qubit (S) coupled on a non-Markovian scale 
to a two-level fluctuator (E), which in turn interacts in 
a Markovian way with a bosonic bath (B) has been de- 
scribed in detail in 1281. 
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Henceforth, for describing the method we will drop the 
subscript to the quantum map F(t) and tacitly assume 
we refer to the smallest embedding such that the map is 
Markovian and governed by Eqn. ([9]). 

Moreover, if the Hamiltonian is composed of the drift 
term Hd and control terms Hj with piecewise constant 
control amplitudes Uj(t k ) for tk G [0,T] 

H{t k ) := Hd + ^UjitkjHj with Uj{t k ) £«CM (10) 
o 

then Eqn. ([9]) defines a bilinear control system. 

With these stipulations, the GRAPE algorithm can be 
lifted to the superoperator level in order to cope with 
open systems by numerically optimising the trace fidelity 

/ tr := Re tr{Ad^ target F(T)} (11) 

for fixed final time T. For simplicity, we henceforth as- 
sume equal time spacing At := t k — t k -\ for all time 
slots k = 1,2,..., M, so T = M ■ At. Therefore 
F(T) = F M ■ F M -i ■■■F k ---F 2 -F 1 with every map tak- 
ing the form F k = exp{ — (i a,d H n k ) +T(t k ))At} leads to 
the derivatives 

4^ = - Re tr { Adl -F M ■ F M ^ ■ ■ ■ F k+1 x 

x(i^ Hj + ^^l)F k AtxF k . 1 ...F 2 .F 1 } 

(12) 

for the recursive gradient scheme 

^ H " 1) (**)=«f(t*)+«r4fe , (13) 

where often the uniform At is absorbed into the step size 
a r > 0. It gives the update from iteration r to r + 1 of 
the control amplitude Uj to control Hj in time slot t k . 

Numerical Setting 

Numerical openGRAPE typically started from some 50 
initial conditions to each fixed final time taking then some 
r = 10 — 30 x 10 3 iterations (see Eqn. [13]) to arrive at one 
point in the top curve shown as upper trace in Fig. [3] 

In contrast, for finding time-optimised controls in the 
closed reference system, we used GRAPE for tracking top 
curves: this is done by performing optimisations with 
fixed final time, which is then successively decreased so 
as to give a top curve g{T) of quality against duration 
of control, a standard procedure used in, e.g., Ref. [12] • 
Finding controls for each fixed final time was typically 
starting out from some 20 random initial control se- 
quences. Convergence to one of the points in time (where 
Fig. [3] shows mean and extremes for a familiy of 15 dif- 
ferent such optimised control sequences) required some 
r = 1000 recursive iterations each. — Numerical exper- 
iments were carried out on single workstations with 512 
MHz to 1.2 GHz tact rates and 512 MB RAM. 

Clearly, there is no guarantee of finding the global op- 
timum this way, yet the improvements are substantial. 



III. EXPLORING APPLICATIONS BY MODEL 
SYSTEMS 

By way of example, the purpose of this section is to 
demonstrate the power of optimal control of open quan- 
tum systems as a realistic means for protecting from re- 
laxation. In order to compare the results with idealised 
scenarios of 'decoherence-free subspaces' and 'bang-bang 
decoupling', we choose two model systems that can par- 
tially be tracted by algebraic means. Comparing nu- 
merical results with analytical ones will thus elucidate 
the pros of numerical optimal control over previous ap- 
proaches. — In order to avoid misunderstandings, how- 
ever, we should emphasize our algorithmic approach to 
controlling open systems (openGRAPE) is by no means 
limited to operating within such predesigned subspaces 
of weak decoherence: e.g., in Ref. [28| we have worked in 
the full Liouville space of a non-Markovian target system. 
Yet, not only are subspaces of weak decoherence prac- 
tically important, they also lend themselves to demon- 
strate the advantages of relaxation-optimised control in 
the case of Markovian systems with time independent 
relaxation operator T, which we focus on in this section. 

The starting point is the usual encoding of one logical 
qubit in Bell states of two physical ones 

\0)l ^7f{|01) + |10)} = |^ + ) 

|l>L:=7i{|01>-|10}} = hn 

Four elements then span a Hermitian operator subspace 
protected against X" 2 -type relaxation 

E-spanjaiV^XV^I} • (15) 

This can readily be seen, since for any p G B 

T Q (p):=[zz,[zz,p]] = Q , (16) 

where henceforth we use the short-hand zz := er 2 <g> a z /2 
and likewise XX cLS well as lfxi/1 := j I2 ® 0|, ® <7„ <g) 1 2 
for /i, v G {x,y,z,t}. Interpreting Eqn. [16] as perfect 
protection against T 2 -type decoherence is in line with the 
slow-tumbling limit of the Bloch-Redfield relaxation by 
the spin tensor A 2j (o,o) : = 75 (§ zz ~ W2) 

!> 2 (p) := [A\ (q q) , [A 2 (0) o) > P]] = M I zz ' t zz ' P\\ = • 

(17) 

For the sake of being more realistic, the model relax- 
ation superoperator mimicking dipole-dipole relaxation 
within the two spin pairs in the sense of Bloch-Redfield 
theory is extended from covering solely T 2 -type decoher- 
ence to mildly includin g T\ dissipation by taking (for each 
basis state p) the sum [44j 
1 

T(p) := J2 [ A l( mum2 )> [\( mitma )>P]] . ( 18 ) 

mi,m2- — 1 

in which the zeroth-order tensor A 2 ,(o,o) ~ zz ls then 
scaled 100 times stronger than the new terms. So the 
resulting model relaxation rate constants finally become 
T.2 1 : Tf 1 = 4.027 s" 1 : 0.024 s" 1 ~ 170 : 1. 
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A. Controllability Combined with Protectability 
against Relaxation 

In practical applications of a given system, a central 
problem boils down to simultaneously solving two ques- 
tions: (i) is the (sub) system fully controllable and (ii) can 
the (sub) system be decoupled from fast relaxing modes 
while being steered to the target. 

It is for answering these questions in algebraic terms 
that we have chosen the following coupling interactions: 
if the two physical qubits are coupled by a Heisenberg-XX 
interaction and the controls take the form of z-pulses 
acting jointly on the two qubits with opposite sign, one 
obtains the usual fully controllable logical single qubit 
over B, because 

((zl-lz),(xx + yy)) Ue ^ p su(2) , (19) 

where (-) Lie denotes the Lie closure under commutation 
(which here gives (yx — xy) as third generator to su(2)). 

Model System I 
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a,d H h-> ( ad H ) r 

Figure 2: (Colour online) In a physical four-qubit system 
for encoding two logical qubits, the Hamiltonians (in their 
superoperator representations of ad.tr) take the form of 
256 x 256 matrices. In the eigenbasis of T (Eqn. [I8j, the 
drift Hamiltonian adH D1 of System / block diagonalises into 
slowly relaxing modes (blue) with relaxation rate constants 
in the interval [0 s" 1 , 0.060 s~ ], moderately relaxing modes 
(magenta) with [4.01 s _1 , 4.06 s _1 ], and fast relaxing modes 
(red) with [8.02 s" 1 , 8.06 s -1 ]. In System lithe Hamiltonian 
adjf D1+D2 comprises off-diagonal blocks (empty boxes) 
that make the protected modes exchange with the fast 
decaying ones. [NB: for pure T2-relaxation (Eqn. 1 1 T jl . the 
relaxation-rate eigenvalues would further degenerate to 0s _1 
('decoherence-free'), 4s _1 (medium) and 8s _1 (fast) while 
maintaining the same block structure]. 



By coupling two of the above qubit pairs with an Ising- 
ZZ interaction as in Refs. [H, Ho, 53] one gets the stan- 
dard logical two-spin system serving as our reference Sys- 
tem I. it is defined by the drift Hamiltonian Hdi and the 
control Hamiltonians Hci , Hc2 

H D i -=Jxx {xxlt + tlxx + yytl + \lyy) + J zz 1zz\ 
H C i :=zlll - Izll 

H C 2 :=ttzt - tttz , (20) 

where the coupling constants are set to J xx = 2 Hz and 
J zz = 1 Hz. Hence, over the T 2 -decoherence protected 
subspace spanned by the four-qubit Bell basis B®B one 
obtains a fully controllable logical two-qubit system 

(H m ,HcuH C 2)^\ mB = P 5U(4) . (21) 

As illustrated in Fig. [21 in the eigenbasis of T (of Eqn. [18]) 
the Hamiltonian superoperators ad# take block diagonal 
form, where the first block acts on the Liouville subspace 
B®B spanning the states protected against T2-type re- 
laxation. Thus in more abstract terms (and recalling 
Eqn. [3l), the Hamiltonians of System I restricted to the 
T 2 -protected block, {adff D1 , ad^ D1 , ad ffm }\ mB , gener- 
ate Adg(7(4) as group of inner automorphisms over the 
protected states. 



Model System II 

Now, by extending the Ising-ZZ coupling between the 
two qubit pairs to an isotropic Heisenberg-XXX interac- 
tion, one gets what we define as System II. Its drift term 



with the coupling constants being set to J xx — 2 Hz and 
Jxyz = 1 Hz reads 

Hdi+D2 ■= Jxx (xxlt + llxx + yylt + llyy) 
+J xyz (Ixxt + tyyl + Izzl) 

and it takes the system out of the decoherence-protected 
subspace due to the off-diagonal blocks in Fig. EJ so the 
dynamics finds its Lie closure in a much larger algebra 
isomorphic to so(12), 

dim ({H D1+D2 ), Hci,H C 2) Lic = 66 , (23) 

to which su(4) is but a subalgebra. 

Note that e~ t7THc " {H D1+D2 )e™ Hc " = H D1 - D2 for ei- 
ther ^ = 1,2. So invoking Trotter's formula 

lim ( e -*( H Dl + 02)/{2n) e -i(H D l-n2)/(2n)\ri = ^iHoi 
n — >oo 

(24) 

it is easy to see that the dynamics of System II may re- 
duce to the subspace of System I in the limit of infinitely 
many switchings of controls Hci or Hc2 and free evo- 
lution under Hdi + d2- It is in this decoupling limit that 
System II encodes a fully controllable logical two-qubit 
system over the then dynamically protectable basis states 
of B ® B. 

In the following paragraph we may thus compare the 
numerical results of decoherence-protection by optimal 
control with alternative pulse sequences derived by pa- 
per and pen exploiting the Trotter limit. As an example 
we choose the CNOT gate in a logical two-qubit system 
encoded in the protected four-qubit physical basis B®B. 
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Figure 3: (Colour online) Fidelity of a CNOT gate encoded in an open system of four physical qubits in different scenarios of 
System II (see text). For reference, the top panel of (a) shows the top curve go{T) consisting of maximum obtainable fidelities 
against fixed final time T in the absence of relaxation; mean and rmsd are shown for families of 15 independent control sequences 
generated for each T. The lower panel of (a) shows their performance vastly scattering in the presence of relaxation (•) with the 
intervals giving mean ± rmsd for all the 15 control sequences tested (dots for best and worst values), while numerical optimal 
control under explicit relaxation (•) is far superior. In (b) these results are compared with (o) by naive Trotter calculations 
assuming to every interaction the inverse is directly obtainable; (•) depicts a realistic Trotter approach, where the inverse has to 
be explicitly generated. The numbers in brackets (ni , H2,p) give the expansion coefficients n\ , ri2 of Fig. [5] and the max. control 
power p (counted as number of 27r-rotations per second). Note that numerical optimal control requires some four orders of 
magnitude less power than the bang-bang type decoupling from fast relaxing modes used in the Trotter expansions. The upper 
quality limit is imposed by slow Ti-type relaxation (see text). Without relaxation, all the Trotter sequences would achieve 
fidelities between 93 and 99 %, except (*) the ones limited to control fields of powers p < 50 Hz: they would fall below 5%. 



B. Results on Performing Target Operations under 
Simultaneous Decoupling 

The model systems are completely parameterised by 
their respective Master equations, i.e. by putting to- 
gether the Hamiltonian parts of Eqns. ([20)1 for System I 
or Eqn. l(22|) for System II and the relaxative part ex- 
pressed in Eqn. (fl8)l . We will thus compare different 
scenarios of approximating the logical CNOT target gate 
(Ad[/ CNOT ) by the respective quantum map F(T) while 
at the same time, the logical two-qubit subsystem has to 
be decoupled from the fast decaying modes in order to 
remain within a weakly relaxing subspace. This is what 
makes it a demanding simultaneous optimisation task. 
— The numerical and analytical results are summerised 
in Fig. [3j they come about as follows. 



1. Comparison of Relaxation- Optimised and Near 
Time- Optimal Controls 



as bad as giving 65%. With openGRAPE performing 
about two standard deviations better than the mean ob- 
tained without taking relaxation into account, only 2.5% 
of near time-optimal control sequences would roughly be 
expected to reach a fidelity beyond 95% just by chance. 
Fig. [4] then elucidates how the new decoherence avoid- 
ing controls keep the system almost perfectly within 
the slowly-relaxing subspace, whereas conventional near 
time-optimal controls partly sweep through the fast- 
relaxing subspace thus leading to inferior quality. 



(a) (b) 




With decoherence-avoiding numerically optimised con- 
trols one obtains a fidelity beyond 95%, while near time- 
optimal controls show a broad scattering as soon as re- 
laxation is taken into account: among the family of 15 
sequences generated, serendipity may help some of them 
to reach a quality of 85 to 90%, while others perform 



Figure 4: (Colour online) (a) Time evolution of all the pro- 
tected basis states under a typical time-optimised control of 
FigOD Projections into the slowly-relaxing and fast-relaxing 
parts of the Liouville space are shown, (b) Same for the new 
decoherence-avoiding controls. System II (see text) then stays 
almost entirely within the T2-protected subspace. 



(a) (b) (c) (d) 




Figure 5: Impractical academic alternative to numerical control: the algebraic derivation of controls for a CNOT in the slowly- 
relaxing subspace proceeds from (a) the logical two-qubit system via (b) the encoded schematic physical four-qubit system to 
the physical realisations in the settings of (c) System I and (d) System II as defined in the text. Note that by the definition of 
the encoding Bell states (Eqn. [14]) all angles halve upon going from (a) to (b). Effective Hamiltonians tH are represented by 
large frames; ni,»i2 indicate repetitions for the respective Trotter expansions. Black bars are local f-pulses with phases given 
as subscript unless other rotation angles given on top; empty bars denote local z-rotations by flip angle n. These pulses can be 
seen as bang-bang type controls; due to the high repetition rates (ni = 2 but n 2 = 64) the Trotter expansions accumulate field 
strengths of p ~ 500 kHz, while the pulse shapes from optimal control require p < 50 Hz. Expanding the interior Hamiltonian 
in (b) is even more complicated (not shown here) 



2. Comparison to Paper-and-Pen Solutions 

Algebraic alternatives to numerical methods of optimal 
control exploit Trotter's formula for remaining within the 
slowly-relaxing subspace when realising the target, see, 
e.g. [48j]. Though straightforward, they soon become 
unhandy as shown in Fig. [U Assuming for the moment 
that to any evolution under a drift Hd the inverse evo- 
lution under — Hd is directly available, the correspond- 
ing "naive" expansions take almost 3 times the length of 
the numerical results, yet requiring much stronger control 
fields (1 — 17 kHz instead of 50 Hz) as shown in Fig. [3j 
In practice, however, the inverse is often not immedi- 
ately reachable, but will require waiting for periodicity. 
For instance, in the Trotter decomposition of Fig. 03(c), 
the Ising term Hzz '■= tzzl as part of the drift Hamil- 
tonian Hdi+d2 is also needed with negative sign so that 
all terms governed by J xyz in Eqn. [22] cancel and only 
the Heisenberg-XX terms governed by Jxx survive. But 
Hzz cannot be sign-reversed directly by the z-controls 
in the sense tzzl i— ► —tzzl since it clearly commutes 
with the z-controls. Thus one will have to choose evolu- 
tion times (t2 in Fig. [5]) long enough to exploit (quasi) 
periodicity. However, Hdi+D2 shows eigenvalues lacking 
periodicity within practical ranges altogether. Moreover, 
the non-zero eigenvalues of -£/di+,D2 do not even occur in 
pairs of opposite sign, hence there is no unitary transform 
U : H h- > —H = UAW to reverse them, and a forteriori 
there is no local control that could do so either [4fi| ]. 

Yet, when shifting the coupling to J xx = 2.23 Hz to 



introduce a favourable quasi-periodicity, one obtains al- 
most perfect projection (/ tr > 1 — 1CP 10 ) onto the inverse 

drift evolution of System II, to wit U' 1 := e +t i HDl+D2 
after 3.98 sec and onto — U^ 1 after 1.99 sec. Thus the 
identity Ad(_j/-i) = Adm-n may be exploited to cut 
the duration for implementing Ad^-i to 1.99 sec. Yet, 
even with these facilitations, the total length required for 
a realistic Trotter decomposition (with an overall trace 
fidelity of /tr > 94.1 % in the absence of decoherence) 
amounts to some 28.5 sec as shown in Fig.[3j Moreover, 
as soon as one includes very mild Ti-type processes, the 
relaxation rate constants in the decoherence-protected 
subspace are no longer strictly zero (as for pure T2-type 
relaxation), but cover the interval [0s _1 , 0.060s -1 ]. Un- 
der these realistic conditions, a Trotter expansion gives 
no more than 15% fidelity, while the new numerical meth- 
ods allow for realisations beyond 95% fidelity in the same 
setting (even with the original parameter J xx = 2.0 Hz). 



IV. DISCUSSION 

In order to extract strategies of how to fight relaxation 
by means of optimal control, we classify open quantum 
systems (i) by their dynamics being Markovian or non- 
Markovian and (ii) by the (Liouville) state space directly 
representing logical qubits either directly without encod- 
ing or indirectly with one logical qubit being encoded by 
several physical ones. So the subsequent discussion will 
lead to assigning different potential gains to different see- 
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narios as summerised in Table [H 

Before going into them in more detail, recall (from the 
section on numerical setting) that the top curve go(T) 
shall denote the maximum fidelity against final times T 
as obtained for the analogous closed quantum system 
(i.e. setting V = 0) by way of numerical optimal con- 
trol. Moreover, define T* as the smallest time such that 
9o(T*) = 1 — e, where e denotes some error-correction 
threshold. 

First (I) , consider the simple case of a Markovian quan- 
tum system with no encoding between logical and phys- 
ical qubits, and assume go(T) has already been deter- 
mined. If as a trivial instance (La) one had a uniform 
decay rate constant 7 so T = 7I, then the fidelity in 
the presence of relaxation would simply boil down to 
f(T) = g (T) ■ e~ T -\ Define := argmax{/(T)} and 
pick the set of controls leading to go calculated in 
the absence of relaxation for tracking go (T) . In the sim- 
plest setting, they would already be 'optimal' without 
ever having resorted to optimising an explicitly open sys- 
tem. More roughly, the time-optimal controls at T = T„ 
already provide a good approximation to fighting relax- 
ation if T* — > is small, i.e. if 7 > is small. 

Next consider a Markovian system without coding, 
where T 7^ 7I is not fully degenerate (Lb). Let {7-,} 
denote the set of (the real parts of the) eigenvalues of 
r. Then, by convexity of {e -7 * | t, 7 > 0}, the following 
rough [53j yet useful limits to the fidelity f(T) obtainable 
in the open system apply 

T N 2 N 2 

<?o(T)-ex P {-— J» < f(T) < go (T)- — Y,^ T ■ 

i=l i=i 

(25) 

Hence the optimisation task in the open system 
amounts to approximating the target unitary gate 
(Adf7 targot ) by the quantum map F(T) resulting from evo- 
lution under the controls subject to the condition that 
modes of different decay rate constants 7? ^ -fk are in- 
terchanged to the least possible amount during the en- 
tire duration < t < T of the controls. An application of 
this strategy known in NMR spectroscopy as TROSY [H3| 
makes use of differential line broadening [5l| and partial 
cancellation of relaxative contributions. Clearly, unless 
the eigenvalues 7,' do not significantly disperse, the ad- 
vantage by optimal control under explicit relaxation will 
be modest, since the potential gain in this scenario relates 
to the variance a ({jj}). 

The situation becomes significantly more rewarding 
when moving to the category (II) of optimisations re- 
stricted to a weakly relaxing (physical) subspace used 
to encode logical qubits. A focus of this work has been 
on showing that for Markovian systems encoding logical 
qubits, the knowledge of the relaxation parameters trans- 
lates into significant advantages of relaxation-optimised 
controls over time-optimised ones. This is due to a dual 
effect: openGRAPE readily decouples the encoding sub- 
system from fast relaxing modes while simultaneously 
generating a quantum map of (close to) best match to 



Table I: Gain Potential for Relaxation-Optimised Controls 
versus Time-Optimised Controls 



Category 


Markovian 


non-Mar kovian 


encoding: 






protected subspace 


big 


(difficult") 


no encoding: 






full Liouville space 


small-medium 


medium-big [28] 



"The problem actually roots in finding a viable protected subspace 
rather than drawing profit from it. 

the target unitary. Clearly, the more the decay of the 
subspace differs from its embedding, the larger the ad- 
vantage of relaxation-optimised control becomes. More- 
over, as soon as the relaxation-rate constants of the pro- 
tected subsystem also disperse among themselves, modes 
of different decay should again only be interchanged to 
the least amount necessary — thus elucidating the very in- 
tricate interplay of simultaneous optimisation tasks that 
makes them prone for numerical strategies. 

In contrast, in the case of entirely unknown relax- 
ation characteristics, where, e.g., model building and sys- 
tem identification of the relaxative part is precluded or 
too costly, we have demonstrated that guesses of time- 
optimal control sequences as obtained from the analo- 
gous closed system may — just by chance — cope with re- 
laxation. This comes at the cost of making sure a suffi- 
ciently large family of time-optimal controls is ultimately 
tested in the actual experiment for selecting among many 
such candidates by trial and error — clearly no more than 
the second best choice after optimal control under explic- 
itly known relaxation. 

In the non-Markovian case, however, it becomes in gen- 
eral very difficult to find a common weakly relaxing sub- 
space for encoding (II. b): there is no Master equation 
of GKS-Lindblad form, the T(t) of which could serve as 
a guideline to finding protected subspaces. Rather, one 
would have to analyse the corresponding non-Markovian 
Kraus maps for weakly contracted subspaces allowing for 
encodings. — However, in non-Markovian scenarios, the 
pros of relaxation-optimised control already become sig- 
nificant without encoding as has been demonstrated in 



Simultaneous Transfer in Spectroscopy 

Finally, note that the presented algorithm also solves 
(as a by-product) the problem of simultaneous state-to- 
state transfer that may be of interest in coherent spec- 
troscopy While Eqns. l9l and fTTI refer to the full-rank 
linear image F, one may readily project onto the states 
of concern by the appropriate projector II to obtain the 
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respective dynamics and quality factor of the subsystem 
HF = -n(iad H +r) o F (26) 

/tr n) = JOT Re tr { nt ^axgetn F(T) } (27) 

reproducing Eqn.[8]in the limit of II being a rank-1 pro- 
jector. While such rank-1 problems under relaxation were 
treated in the algorithmic setting of openGRAPE put 
forward here allows for projectors of arbitrary rank, e.g., 
1 < rkll < N for n spin-^ qubits with TV := 2™. Clearly, 
the rank equals the number of orthogonal state-to-state 
optimisation problems to be solved simultaneously. 

V. CONCLUSIONS AND OUTLOOK 

We have provided numerical optimal-control tools to 
systematically find near optimal approximations to uni- 
tary target modules in open quantum systems. The pros 
of relaxation-optimised controls over time-optimised ones 
depend on the specific experimental scenario. We have 
extensively discussed strategies for fighting relaxation in 
Markovian and non-Markovian settings with and without 
encoding logical qubits in protected subspaces. Numeri- 
cal results have been complemented by algebraic analysis 
of controllability in protected subspaces under simulta- 
neous decoupling from fast relaxing modes. 

To complement the account on non-Markovian systems 
in [1^1, the progress is quantitatively exemplified in a 
typical Markovian model system of four physical qubits 
encoding two logical ones: when the Master equation is 
known, the new method is systematic and significantly 
superior to near time-optimal realisations, which in turn 
are but a guess when the relaxation process cannot be 
quantitatively characterised. In this case, testing a set of 
10 — 20 such near time-optimal control sequences empiri- 
cally is required for getting acceptable results with more 
confidence, yet on the basis of trial and error. As follows 
by controllability analysis, Trotter-type expansions allow 
for realisations within slowly-relaxing subspaces in the 



limit of infinitely many switchings. However, in realistic 
settings for obtaining inverse interactions, they become 
so lengthy that they only work in the idealised limit of 
both T2 and Ti-decoherence-free subspaces, but fail as 
soon as very mild Ti-relaxation processes occur. 

Optimal control tools like openGRAPE are therefore the 
method of choice in systems with known relaxation pa- 
rameters. They accomplish decoupling from fast relaxing 
modes with several orders of magnitude less decoupling 
power than by typical bang-bang controls. Being appli- 
cable to spin and pseudo-spin systems, they are antici- 
pated to find broad use for fighting relaxation in practical 
quantum control. In a wide range of settings the benefit 
is most prominent when encoding the logical system in a 
protected subspace of a larger physical system. However, 
the situation changes upon shifting to a timevarying T(t) 
[2S| |. or to more advanced non-Markovian models with 
r(u(t)) depending on time via the control amplitudes 
u(t) on timescales comparable to the quantum dynamic 
process. Then the pros of optimal control extend to the 
entire Liouville space, as shown in [2g| . 

In order to fully exploit the power of optimal control 
of open systems the challenge is shifted to (i) thoroughly 
understanding the relaxation mechanisms pertinent to a 
concrete quantum hardware architecture and (ii) being 
able to determine its relaxation parameters to sufficient 
accuracy. 
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